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Abstract 

The results of a search for charged Higgs bosons are presented. The analysis is based on 
4.6 fb" 1 of proton-proton collision data at y/s = 7 TeV collected by the ATLAS experiment at the Large 
Hadron Collider, using top quark pair events with a r lepton in the final state. The data are consistent 
with the expected background from Standard Model processes. Assuming that the branching ratio 
of the charged Higgs boson to a r lepton and a neutrino is 100%, this leads to upper limits on the 
branching ratio of top quark decays to a b quark and a charged Higgs boson between 5% and 1% for 
charged Higgs boson masses ranging from 90 GeV to 160 GeV, respectively. In the context of the 
m^ ax scenario of the MSSM, tan/3 above 12-26, as well as between 1 and 2-6, can be excluded for 
charged Higgs boson masses between 90 GeV and 150 GeV. 
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Abstract: The results of a search for charged Higgs bosons are presented. The analysis 
is based on 4.6 fb" 1 of proton-proton collision data at \fs = 7 TeV collected by the ATLAS 
experiment at the Large Hadron Collider, using top quark pair events with a r lepton in the 
final state. The data are consistent with the expected background from Standard Model 
processes. Assuming that the branching ratio of the charged Higgs boson to a r lepton and 
a neutrino is 100%, this leads to upper limits on the branching ratio of top quark decays to 
a b quark and a charged Higgs boson between 5% and 1% for charged Higgs boson masses 
ranging from 90 GeV to 160 GeV, respectively. In the context of the m™ ax scenario of 
the MSSM, tan/3 above 12-26, as well as between 1 and 2-6, can be excluded for charged 
Higgs boson masses between 90 GeV and 150 GeV. 
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1 Introduction 



Charged Higgs bosons (H + , H~) are predicted by several non-minimal Higgs scenarios, 
such as Two Higgs Doublet Models (2HDM) [1] or models containing Higgs triplets [2-6]. 
As the Standard Model (SM) does not contain any elementary charged scalar particle, the 
observation of a charged Higgs boson 1 would clearly indicate new physics beyond the SM. 
For instance, supersymmetric models predict the existence of charged Higgs bosons. In 
a type-II 2HDM, such as the Higgs sector of the Minimal Supersymmetric extension of 
the Standard Model (MSSM) [7-11], the main H + production mode at the Large Hadron 
Collider (LHC) is through top quark decays t — > bH + , for charged Higgs boson masses 
(m H +) smaller than the top quark mass (m t0 p)- The dominant source of top quarks at 
the LHC is through ti production. The cross section for H + production from single top 
quark events is much smaller and is not considered here. For tan/3 > 2, where tan/3 is the 
ratio of the vacuum expectation values of the two Higgs doublets, the charged Higgs boson 
decay via — > tv is dominant and remains sizeable for 1 < tan/3 < 2 [12]. In this pa- 
per, B(H + — > tv) = 100% is assumed, unless otherwise specified. Under this assumption, 
the combined LEP lower limit for the charged Higgs boson mass is about 90 GeV [13]. 
The Tevatron experiments placed upper limits on B(t — > bH + ) in the 15-20% range for 
m n + < m top [14, 15]. 

This paper describes a search for charged Higgs bosons with masses in the range 90- 
160 GeV, using ti events with a leptonically or hadronically decaying r lepton in the final 
state, i.e. with the topology shown in Fig. 1. Charged Higgs bosons are searched for in a 
model-independent way, hence exclusion limits are given in terms of B(t — > bH + ), as well 
as in the m™ ax scenario [16] of the MSSM. The results are based on 4.6 fb _1 of data from 
pp collisions at y/s = 7 TeV, collected in 2011 with the ATLAS experiment [17] at the 
LHC. Three final states, which are expected to yield the highest sensitivity, are analysed: 

• lepton+jets: ti — > bbWH + — > bb(qq')(T\ ep v), i.e. W decays hadronically and r decays 
into an electron or a muon, with two neutrinos; 

• r+lepton: ti — > bbWH + — > bb(lv)(Tha,dv), i.e. W decays leptonically (with I = e, [/,) 
and r decays hadronically; 

• r+jets: ti — > bbWH + — > bb{qq'){T^.^v), i.e. both W and r decay hadronically. 

In Section 2, the data and simulated samples used in this analysis are described. In 
Section 3, the reconstruction of physics objects in ATLAS is discussed. Sections 4-6 present 
results obtained in the lepton+jets, r+lepton and r+jets channels, respectively. Systematic 
uncertainties are discussed in Section 7, before exclusion limits in terms of B(t — > bH + ) 
and tan /3 are presented in Section 8. Finally, a summary is given in Section 9. 

1 In the following, charged Higgs bosons are denoted H + , with the charge-conjugate H~ always implied. 
Hence, r denotes a positively charged r lepton. 
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Figure 1. Example of a leading-order Feynman diagram for the production of tt events arising from 
gluon fusion, where a top quark decays to a charged Higgs boson, followed by the decay H + — > tv. 

2 Data and simulated events 

The ATLAS detector [17] consists of an inner tracking detector with a coverage in pseudo- 
rapidity 2 up to \r]\ = 2.5, surrounded by a thin 2 T superconducting solenoid, a calorimeter 
system extending upto \rj\ = 4.9 for the detection of electrons, photons and hadronic jets, 
and a large muon spectrometer extending up to |?7| = 2.7 that measures the deflection of 
muon tracks in the field of three superconducting toroid magnets. A three-level trigger sys- 
tem is used. The first level trigger is implemented in hardware, using a subset of detector 
information to reduce the event rate to a design value of at most 75 kHz. This is followed 
by two software-based trigger levels, which together reduce the event rate to about 300 Hz. 

Only data taken with all ATLAS sub-systems operational are used. It results in an in- 
tegrated luminosity of 4.6 fb _1 for the 2011 data-taking period. The integrated luminosity 
has an uncertainty of 3.9%, measured as described in Refs. [18, 19] and based on the whole 
2011 dataset. Following basic data quality checks, further event cleaning is performed by 
demanding that no jet is consistent with having originated from instrumental effects, such 
as large noise signals in one or several channels of the hadronic end-cap calorimeter, co- 
herent noise in the electromagnetic calorimeter, or non-collision backgrounds. In addition, 
events are discarded if the reconstructed vertex with the largest sum of squared track mo- 
menta has fewer than five associated tracks with transverse momenta > 400 MeV. 

The background processes that enter this search include the SM pair production of top 
quarks ti — > bbW + W~ , as well as the production of single top quark, VF+jets, Z/7*+jets, 
diboson and multi-jet events. Data-driven methods are used in order to estimate the multi- 
jet background, as well as the backgrounds with intrinsic missing transverse momentum 
where electrons or jets are misidentified as hadronically decaying r leptons. The modelling 

2 ATLAS uses a right-handed coordinate system with its origin at the nominal interaction point (IP) in 
the centre of the detector and the z-axis along the beam pipe. The a:-axis points from the IP to the centre 
of the LHC ring, and the y-axis points upwards. Cylindrical coordinates (r, (j>) are used in the transverse 
plane, <f> being the azimuthal angle around the beam pipe. The pseudorapidity is defined in terms of the 
polar angle 8 as rj = — lntan(#/2). 
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of SM it and single top quark events is performed with MC@NLO [20], except for the 
t-channel single top quark production where AcerMC [21] is used. The top quark mass is 
set to 172.5 GeV and the set of parton distribution functions used is CT10 [22]. For the 
events generated with MC@NLO, the parton shower, hadronisation and underlying event 
are added using HERWIG [23] and JIMMY [24]. PYTHIA [25] is instead used for events 
generated with AcerMC. Inclusive cross sections are taken from the approximate next- 
to-next-to-leading-order (NNLO) predictions for ti production [26], for single top quark 
production in the t-channel and s-channel [27, 28], as well as for Wt production [29]. Over- 
laps between Wt and SM tt final states are removed [30]. Single vector boson (W and 
Z/Y) production is simulated with ALPGEN [31] interfaced to HERWIG and JIMMY, 
using CTEQ6.1 [32] parton distribution functions. The additional partons produced in the 
matrix element part of the event generation can be light partons or heavy quarks. In the 
latter case, dedicated samples with matrix elements for the production of massive bb or cc 
pairs are used. Diboson events (VFVK, WZ and ZZ) are generated using HERWIG. The 
cross sections are normalised to NNLO predictions for W and Z/7* production [33, 34] 
and to next-to- leading-order (NLO) predictions for diboson production [35]. 

The SM background samples are summarised in Table 1. In addition, three types of 
signal samples are produced with PYTHIA for 90 GeV < m H + < 160 GeV: tt — > bbH + W~, 
ti — > bbH~W + and tt — > bbH + H~, where the charged Higgs bosons decay as H + — > tu. 
The cross section for each of these three processes depends only on the total ti production 
cross section (167 pb) and the branching ratio B(t — > bH + ). TAUOLA [36] is used for r 
decays, and PHOTOS [37] is used for photon radiation from charged leptons. 



Process 


Generator Cross section [pb] 


SM ti with at least one lepton £ = e, /j,, t 


MC@NLO [20] 91 [26] 


Single top quark t-channel (with €) 
Single top quark s-channel (with £) 
Single top quark Wi-channel (inclusive) 


AcerMC [21] 21 [27] 
MC@NLO [20] 1.5 [28] 
MC@NLO [20] 16 [29] 


W^£u 


ALPGEN [31] 3.1 x 10 4 [33] 


Z/Y U with m{U) > 10 GeV 


ALPGEN [31] 1.5 x 10 4 [34] 


WW 

zz 
wz 


HERWIG [23] 17 [35] 
HERWIG [23] 1.3 [35] 
HERWIG [23] 5.5 [35] 


H + signal with B(t bH + ) = 5% 


PYTHIA [25] 16 



Table 1. Cross sections for the simulated processes and generators used to model them. 



The event generators are tuned in order to describe the ATLAS data. The parameter 
sets AUET2 [38] and AUET2B [39] are used for events for which hadronisation is simulated 
using HERWIG/ JIMMY and PYTHIA, respectively. To take into account the presence of 
multiple interactions (around nine, on average) occurring in the same and neighbouring 
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bunch crossings (referred to as pile-up), simulated minimum bias events are added to the 
hard process in each generated event. Prior to the analysis, simulated events are reweighted 
in order to match the distribution of the average number of pile-up interactions in the data. 
All generated events are propagated through a detailed GEANT4 simulation [40, 41] of the 
ATLAS detector and are reconstructed with the same algorithms as the data. 

3 Physics object reconstruction 

3.1 Electrons 

Electrons are reconstructed by matching clustered energy deposits in the electromagnetic 
calorimeter to tracks reconstructed in the inner detector. The electron candidates are 
required to meet quality requirements based on the expected shower shape [42], to have a 
transverse energy > 20 GeV and to be in the fiducial volume of the detector, \rj\ < 2.47 
(the transition region between the barrel and end-cap calorimeters, 1.37 < \rj\ < 1.52, is 
excluded). Additionally, £7t and 77-dependent calorimeter (tracking) isolation requirements 
are imposed in a cone with a radius 3 AR = 0.2 (0.3) around the electron position, excluding 
the electron object itself, with an efficiency of about 90% for true isolated electrons. 

3.2 Muons 

Muon candidates are required to contain matching inner detector and muon spectrometer 
tracks [43], as well as to have pt > 15 GeV and \r]\ < 2.5. Only isolated muons are ac- 
cepted by requiring that the transverse energy deposited in the calorimeters (the transverse 
momentum of the inner detector tracks) in a cone of radius AR = 0.2 (0.3) around the 
muon amounts to less than 4 GeV (2.5 GeV). The energy and momentum of the muon are 
excluded from the cone when applying these isolation requirements. 

3.3 Jets 

Jets are reconstructed using the anti-kt algorithm [44, 45] with a size parameter value of 
R = 0.4. The jet finder uses reconstructed three-dimensional, noise-suppressed clusters of 
calorimeter cells [46] . Jets are calibrated to the hadronic energy scale with correction factors 
based on simulation [47, 48]. A method that allows for the identification and selection of 
jets originating from the hard-scatter interaction through the use of tracking and vertexing 
information is used [49]. This is referred to as the "Jet Vertex Fraction" (JVF), defined as 
the fraction of the total momentum of the charged particle tracks associated to the jet which 
belongs to tracks that are also compatible with the primary vertex. By convention, jets 
with no associated tracks are assigned a JVF value of —1 in order to keep a high efficiency 
for jets at large values of 77, outside the range of the inner tracking detectors. The jet 
selection based on this discriminant is shown to be insensitive to pile-up. A requirement 
of I JVF I > 0.75 is placed on all jets during event selection. In order to identify the jets 
initiated by b quarks, an algorithm is used that combines impact-parameter information 

3 Ai? = \J (A7/) 2 + (A(f>) 2 , where A77 is the difference in pseudorapidity of the two objects in question, 
and A<f> is the difference between their azimuthal angles. 
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with the explicit determination of a secondary vertex [50] . A working point is chosen that 
corresponds to an average efficiency of about 70% for b jets with pr > 20 GeV in tt events 
and a light-quark jet rejection factor of about 130. Since the 6-tagger relies on the inner 
tracking detectors, the acceptance region for jets is restricted to \rj\ < 2.4. 

3.4 r jets 

In order to reconstruct hadronically decaying r leptons, anti-fcj jets with either one or three 
associated tracks reconstructed in the inner detector and depositing Er > 10 GeV in the 
calorimeter are considered as r candidates [51]. Dedicated algorithms are used in order to 
reject electrons and muons. Hadronic r decays are identified using a likelihood criterion 
designed to discriminate against quark- and gluon-initiated jets by using the shower shape 
and tracking variables as inputs. A working point with an efficiency of about 30% for 
hadronically decaying r leptons with px > 20 GeV in Z — >■ tt events is chosen, leading to 
a rejection factor of about 100-1000 for jets. The rejection factor depends on the px and r\ 
of the candidate and the number of associated tracks. The r candidates are further required 
to have a visible transverse momentum of at least 20 GeV and to be within \r\\ < 2.3. The 
selected r candidates are henceforth referred to as "r jets". 

3.5 Removal of geometric overlaps between objects 

When candidates selected using the criteria above overlap geometrically, the following 
procedures are applied, in this order: muon candidates are rejected if they are found 
within AR < 0.4 of any jet with px > 25 GeV; a r jet is rejected if found within AR < 0.2 
of a selected muon or electron; jets are removed if they are within AR < 0.2 of a selected 
r object or electron. 

3.6 Missing transverse momentum 

The missing transverse momentum and its magnitude E™ 1SS [52] are reconstructed from 
three-dimensional, noise-suppressed clusters of cells in the calorimeter and from muon 
tracks reconstructed in the muon spectrometer and the inner tracking detectors. Clusters 
of calorimeter cells belonging to jets (including r jets) with px > 20 GeV are calibrated to 
the hadronic energy scale. Calorimeter cells not associated with any object are also taken 
into account and they are calibrated at the electromagnetic energy scale. In order to deal 
appropriately with the energy deposited by muons in the calorimeters, the contributions 
of muons to E™ 1SS are calculated differently for isolated and non-isolated muons. 

4 Analysis of the lepton+jets channel 

This analysis relies on the detection of lepton+jets decays of tt events, where the charged 
lepton I (electron or muon) arises from H + — > T\ cp v, while the jets arise from a hadronically 
decaying W boson, i.e. tt — > bbWH + — > bb{qq'){T\ ep u). 
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4.1 Event selection 

The lepton+jets analysis uses events passing a single-lepton trigger with an Ey threshold 
of 20-22 GeV for electrons 4 and a p^ threshold of 18 GeV for muons. These thresholds 
are low enough to guarantee that electrons and muons chosen for the analysis are in the 
plateau region of the trigger-efficiency curve. In addition, to select a sample of lepton+jets 
events enriched in ti candidates, the following requirements are applied: 

• exactly one lepton having Et > 25 GeV (electron) or > 20 GeV (muon) and 
matched to the corresponding trigger object, with neither a second lepton nor a r jet 
in the event; 

• at least four jets having p^ > 20 GeV, with exactly two of them being 6-tagged; 

• E™ 1SS > 40 GeV and, in order to discriminate between E™ ss arising from isolated 
neutrinos and from poorly reconstructed leptons, this requirement is tightened to 
E™ 1SS x | sin Ai^missl > 20 GeV if the azimuthal angle A<^ m i SS between the lepton 
and E!f lss is smaller than 7r/6. 

Having selected a lepton+jets sample enriched in ti candidates, jets must be assigned 
correctly to the decay products of each W boson (with a mass my/ = 80.4 GeV) and top 
quark. In particular, the hadronic side of the event is identified by selecting the combination 
of one 6-tagged jet (b) and two untagged jets (j) that minimises: 



where <7t p = 17 GeV and ay/ = 10 GeV are the widths of the reconstructed top quark and 
W boson mass distributions, as measured in simulated ti events. Using information about 
the correctly identified combinations in the generated events, the jet assignment efficiency 
is found to be 72%. Events with x 2 > 5 are rejected in order to select well-reconstructed 
hadronic top quark candidates. 

4.2 Data-driven estimation of backgrounds with misidentified leptons 

While the ATLAS lepton identification gives a very pure sample of candidates, there is a 
non-negligible contribution from non-isolated leptons arising from the semileptonic decay 
of hadrons containing b or c quarks, from the decay-in-flight of tt^ or K mesons and, in 
the case of misidentified electron objects, from the reconstruction of 7r° mesons, photon 
conversions or shower fluctuations. All leptons coming from such mechanisms are referred 
to as misidentified leptons, as opposed to truly isolated leptons (e.g. from the prompt decay 
of W or Z bosons), which are referred to as real leptons. The data-driven estimation of 
the number of misidentified leptons passing the lepton selections of Sections 3.1 and 3.2 is 
based on exploiting differences in the lepton identification between real and misidentified 

4 The electron trigger threshold was increased from 20 GeV to 22 GeV towards the end of data-taking 





in 2011. 



-7- 



electrons or muons. Two data samples are denned, which differ only in the lepton identi- 
fication criteria. The tight sample contains mostly events with real leptons and uses the 
same lepton selection as in the analysis. The loose sample contains mostly events with 
misidentified leptons. This latter sample is obtained by loosening the isolation and identi- 
fication requirements for the leptons. For loose electrons, the isolation requirements have 
an efficiency of about 98% for true isolated electrons, compared to 90% in the tight sample. 
For loose muons, the isolation requirement is removed. By construction, the tight sample 
is therefore a subset of the loose sample. 

Let iV,r and iV^ (iV r T and N^) be the number of events containing real and misidentified 
leptons, respectively, passing a loose (tight) selection. The numbers of events containing 
one loose or tight lepton are given by: 

N L = + iV r L , (4.2) 

N T = N^ + N?. (4.3) 

Defining p T and p m as: 

^ = # ^d p m = ^, (4.4) 

the number of misidentified leptons passing the tight selection can then be written as: 

N T = _PE_ {jhN I._ N T )m (45) 
Pr Pm 

The main ingredients of this data-driven method are thus the relative efficiencies p r 
and p m for a real or a misidentified lepton, respectively, to be detected as a tight lepton. 
The lepton identification efficiency p T is measured using a tag-and-probe method on Z — > 11 
data events with a dilepton invariant mass between 86 GeV and 96 GeV, where one lepton 
is required to fulfill tight selection criteria. The rate at which the other lepton passes 
the same tight selection criteria defines p r . The average values of the electron and muon 
identification efficiencies are 80% and 97%, respectively. On the other hand, a control 
sample with misidentified leptons is selected by considering events in the data with exactly 
one lepton passing the loose criteria. In order to select events dominated by multi-jet 
production, E™ ss is required to be between 5 GeV and 20 GeV. Residual true leptons 
contribute at a level below 10% and are subtracted from this sample using simulation. After 
this subtraction, the rate at which a loose lepton passes tight selection criteria defines the 
misidentification rate p m . The average values of the electron and muon misidentification 
probabilities are 18% and 29%, respectively. In the final parameterisation of p r and p m , 
dependencies on the pseudorapidity of the lepton, its distance AR to the nearest jet and 
the leading jet pt are taken into account. 

4.3 Reconstruction of discriminating variables after the selection cuts 

The analysis uses two variables that discriminate between leptons produced in r — > lv[V T 
and leptons coming directly from W boson decays. The first discriminating variable is the 
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invariant mass iribi of the b jet and the charged lepton I coming from the same top quark 
candidate, or more conveniently, cos Of defined as: 

co< = 2 2ml 2 - 1 * f: P \ - 1- (4-6) 

m top m iy m top m W 

Both 7tt| and are neglected, hence m?, — 2p b -p l , where p b and are the four-momenta 
of the b jet and of the charged lepton I, respectively. The presence of a charged Higgs 
boson in a leptonic top quark decay reduces the invariant product p b ■ p l , when compared 
to VF-mediated top quark decays, leading to cosfy* values closer to —1. 

The second discriminating variable is the transverse mass [53] , obtained by fulfilling 
the constraint (p miss + p l + p b ) 2 = m^ op on the leptonic side of lepton+jets tt events. 
More than one neutrino accounts for the invisible four-momentum p mms and its transverse 
component pr miss - By construction, gives an event-by-event lower bound on the mass 
of the leptonically decaying charged (W or Higgs) boson produced in the top quark decay, 
and it can be written as: 

{m%? = (V<p + (^t +PT b +PT miss ) 2 ~p\j - (pt 1 + PT miss ) 2 ■ (4.7) 

The cos 6i distribution measured in the data is shown in Fig. 2(a), superimposed 
on the predicted background, determined with a data-driven method for the multi-jet 
background and simulation for the other SM backgrounds. In the presence of a charged 
Higgs boson in the top quark decays, with a branching ratio B(t — ?• bH + ), the contribution 
of tt — > bbW + W~ events in the background is scaled according to this branching ratio. A 
control region enriched in tt — > bbW + W~ events is defined by requiring —0.2 < cosfy* < 1. 
In Section 8, this sample is used to fit the branching ratio B(t — > bH + ) and the product 
of the cross section abbww, the luminosity, the selection efficiency and acceptance for 
tt — > bbW + W~ , simultaneously with the likelihood for the signal estimation. In turn, 
this ensures that the final results, and in particular the upper limit on B{t — > bH + ), are 
independent of the assumed theoretical production cross section for tt. With a branching 
fraction B{t — > bH + ) = 5%, the signal contamination in the control region would range 
from 1.3% for m H + = 90 GeV to 0.4% for m H + = 160 GeV. The signal region is defined 
by requiring cos 8* < — 0.6 and rn^ < 60 GeV, where: 

= ^/2^£;- iss (l-cosA^ imiss ). (4.8) 

This is done in order to suppress the background from events with a W boson decaying 
directly into electrons or muons. For events in the signal region, , shown in Fig. 2(b), 
is used as a discriminating variable to search for charged Higgs bosons. Table 2 lists the 
contributions to the signal region of the SM processes and of tt events with at least one decay 
t —> bH + , assuming rrijj+ = 130 GeV and B(t —> bH + ) = 5%. When including signal in the 
prediction, the simulated SM tt contribution is scaled according to this branching ratio. 
The data are consistent with the predicted SM background and no significant deformation 
of the distribution is observed. 
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120 ± 4 ± 25 
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990 ±21 ± 140 



Table 2. Expected event yields in the signal region of the lepton+jets final state, and comparison 
with 4.6 fb" 1 of data. A cross section of 167 pb is assumed for the SM tt background. The numbers 
shown in the last two rows, for a hypothetical H + signal with m#+ = 130 GeV, are obtained with 
B(t — > bH + ) = 5%. Both statistical and systematic uncertainties are shown, in this order. 
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Figure 2. Distribution of (a) cos 8* and (b) m^, in the signal region (cos Of < —0.6, rn^ < 60 GeV) 
for the latter. The dashed line corresponds to the SM-only hypothesis and the hatched area around 
it shows the total uncertainty for the SM backgrounds, where "Others" refers to the contribution 
of all SM processes except tt — > bbW + W~ . The solid line shows the predicted contribution of 
signal+background in the presence of a 130 GeV charged Higgs boson, assuming Bit — > bH + ) = 5% 
and B(H + —> tv) = 100%. The light area below the solid line corresponds to the contribution of 
the H + signal, stacked on top of the scaled tt — > bbW + W~ background and other SM processes. 



5 Analysis of the T+lepton channel 

This analysis relies on the detection of r±lepton decays of tt events, where the hadronically 
decaying r lepton arises from H + — > Thad^i while an electron or muon comes from the decay 
of the W boson, i.e. tt — > bbWH + — > 66(/i/)(r ha d^)- 
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5.1 Event selection 



The r+lepton analysis relies on the same single- lepton trigger signatures as the lepton+jets 
analysis presented in Section 4. In order to select r+lepton events, the following require- 
ments are made: 

• exactly one lepton, having > 25 GeV (electron) or pi > 20 GeV (muon) and 
matched to the corresponding trigger object, and no other electron or muon; 

• exactly one r jet having pr > 20 GeV and an electric charge opposite to that of the 
lepton; 

• at least two jets having p^ > 20 GeV, including at least one 6-tagged jet; 

• Spt > 100 GeV in order to suppress multi-jet events, where YIpt is the sum of the 
transverse momenta of all tracks associated with the primary vertex. Tracks entering 
the sum must pass quality cuts on the number of hits and have pt > 1 GeV. As this 
variable is based on tracks from the primary vertex (as opposed to energy deposits 
in the calorimeter), it is robust against pile- up. 

E™ 1SS is used as the discriminating variable to distinguish between SM tt events and 
those where top quark decays are mediated by a charged Higgs boson, in which case the 
neutrinos are likely to carry away more energy. 

5.2 Data-driven estimation of backgrounds with misidentified leptons 

The estimation of the backgrounds with misidentified leptons uses the data-driven method 
described in Section 4.2. When implementing the method, the dependence of real and 
misidentification rates on the 6-tagged jet multiplicity are taken into account, as well as 
the requirement for one r jet (instead of a r jet veto). 

5.3 Backgrounds with electrons and jets misidentified as r jets 

The background with electrons misidentified as r jets is estimated using a Z — > ee con- 
trol region in the data [51], where one electron is reconstructed as a r jet. The measured 
misidentification probabilities, which have an average value of 0.2%, are then applied to all 
simulated events in the r+lepton analysis. Simulation studies show that this application 
is valid, as the misidentification probabilities for Z — > ee and tt events are similar. 

A data-driven method applied to a control sample enriched in VF+jets events is used 
to measure the probability for a jet to be misidentified as a hadronically decaying r lepton. 
This measured probability is used to predict the yield of background events due to jet— > r 
misidentification. Like jets from the hard process in the dominant it background, jets in the 
control sample originate predominantly from quarks instead of gluons. The main difference 
between ti and VF+jets events is the different fraction of b jets, which is smaller in W+jets 
events. However, the probability for a b jet to be misidentified as a r jet is smaller than 
the corresponding probability for a light-quark jet, because the average track multiplicity 
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is higher for b jets. Moreover, the visible mass measurement used in the r identification 
provides further discrimination between b jets and r jets. Differences in jet composition 
(e.g. the ratio of gluons to quarks) between ti and W+jets, assessed using simulation, are 
taken into account as systematic uncertainties. These also cover the dependence of the 
probability on whether a b jet or a light-quark jet is misidentified as a r jet. Events in 
the control region are required to pass the same single-lepton trigger, data quality and 
lepton requirements as in the r+lepton event selection. Additionally, a r candidate and 
E!f lss > 40 GeV are required, and events with 6-tagged jets are vetoed. Simulated events 
with a true r contribute at a level below 0.5% and are subtracted. The r candidates are 
required to have pi > 20 GeV, |r/| < 2.3, and cannot be within AR = 0.2 of any electron 
or muon. They are also not required to pass r identification. The jet — > r misidentification 
probability is defined as the number of objects passing the full r identification divided by 
the number prior to requiring identification. This misidentification probability is evaluated 
separately for r candidates with one or three associated tracks (the corresponding average 
values are about 7% and 2%, respectively) and, in addition, it is measured as a function of 
both pi and rj. 

In order to predict the background for the charged Higgs boson search, the measured 
jet— > t misidentification probability is applied to simulated ti, single top quark, VF+jets, 
Z/7*+jets and diboson events, all of which are required to pass the full event selection 
except for the r identification. For these events, r candidates not overlapping with a 
true r lepton or a true electron, but otherwise fulfilling the same requirements as in the 
denominator of the misidentification probability, are identified. Each of them is considered 
separately to be potentially misidentified as a r jet. In order to avoid counting the same 
object twice, each jet that corresponds to a r candidate is removed from the event. The 
number of reconstructed jets and the number of 6-tagged jets are adjusted accordingly. If, 
after taking this into consideration, the event passes the r+lepton selection, it is counted as 
a background event with a weight given by the misidentification probability corresponding 
to the pt and r\ of the r candidate. The predicted numbers of events from this data- 
driven method and from simulation are shown in Table 3. The backgrounds arising from 
the jet—)- r misidentification are not well modelled in simulation, which is why they are 
estimated using data-driven methods. 

5.4 Event yields and E™ 1SS distribution after the selection cuts 

Table 4 shows the expected number of background events for the SM-only hypothesis and 
the observation in the data. The total number of predicted events (signal+background) in 
the presence of a 130 GeV charged Higgs boson with B(t — > bH + ) = 5% is also shown. The 
r+lepton analysis relies on the theoretical ti production cross section a t i = 167l}g pb [26] 
for the background estimation. In the presence of a charged Higgs boson in the top quark 
decays, with a branching ratio B(t — > bH + ), the contributions of ti —> bbW + W~ events 
in the backgrounds with true or misidentified r jets are scaled according to this branching 
ratio. The background with correctly reconstructed r jets is obtained with simulation. The 
data are found to be consistent with the expectation for the background-only hypothesis. 
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Sample 


Data-driven method [events] 


Simulation [events] 


tt 


900 ± 15 


877 ±6 


W+jets 


150 ± 3 


145 ±9 


Single top quark 


81 ± 1 


61 ±2 


Z/ 7 *+jets 


44 ± 1 


69 ±4 


Diboson 


6± 1 


8± 1 



Table 3. Application of the misidcntification probability obtained from I^+jets events in the data, 
for the r+lepton channel. The predictions of the background contributions based on data-driven 
misidentification probabilities and on simulation are given, with statistical uncertainties only. In 
both cases, all top quarks are assumed to decay via t —> bW. 

The E™ 1SS distributions for the r + e and r + fi channels, after all selection cuts are applied, 
are shown in Fig. 3. 



Sample 


Event yield (r+lepton) 
r + e r + fx 


True r+lepton 
Misidentified jet— t r 
Misidentified e — > r 
Misidentified leptons 


430 + 14+ 59 570 + 15+ 75 
510 + 23+ 86 660 + 26 + 110 
33+4+ 5 34+4+ 6 
39+ 10+ 20 90 + 10+ 34 


All SM backgrounds 


1010 + 30 + 110 1360 + 30 + 140 


Data 


880 1219 


t -> bH+ (130 GeV) 
Signal+background 


220+ 6+ 29 310+ 7+ 39 
1160 ± 30 ± 100 1570 ± 30 ± 130 



Table 4. Expected event yields after all selection cuts in the r+lepton channel and comparison 
with 4.6 fb" 1 of data. The numbers in the last two rows, obtained for a hypothetical H + signal with 
m H+ = 130 GeV, are obtained with B(t -> bH+) = 5%. All other rows assume B(t -> bW) = 100%. 
Both statistical and systematic uncertainties are shown, in this order. 



6 Analysis of the r+jets channel 

The analysis presented here relies on the detection of r+jets decays of tt events, where 
the hadronically decaying r lepton arises from H + — > t^^v, while the jets come from a 
hadronically decaying W boson, i.e. ti — > bbWH + — > 66(gg')( r had z/ )- 

6.1 Event selection 

The r+jets analysis uses events passing a r + E™ 1SS trigger with a threshold of 29 GeV 
on the r object and 35 GeV on calorimeter-based £™ 1SS . The following requirements are 
applied, in this order: 
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Figure 3. E™ lss distribution after all selection cuts in the r+lepton channel, for (a) r+electron 
and (b) r+muon final states. The dashed line corresponds to the SM-only hypothesis and the 
hatched area around it shows the total uncertainty for the SM backgrounds. The solid line shows 
the predicted contribution of signal+background in the presence of a 130 GeV charged Higgs boson 
with B(t -> bH+) = 5% and B(H+ -> tv) = 100%. The contributions of tt -> bbW+W events in 
the backgrounds with true or misidentified r jets are scaled down accordingly. 



• at least four jets (excluding r jets) having px > 20 GeV, of which at least one is 
6-tagged; 

• exactly one r jet with pT^ > 40 GeV, found within \rj\ < 2.3 and matched to a r 
trigger object; 

• neither a second r jet with p7^ > 20 GeV, nor any electrons with > 20 GeV, nor 
any muons with p^ > 15 GeV; 

• E™ iss > 65 GeV; 

• to reject events in which a large reconstructed i?™ lss is due to the limited resolution of 
the energy measurement, the following ratio based on the ^2pt definition of Section 5 
must satisfy: 

irimiss 

> 13; 



0.5 GeV 1 / 2 • y/Yj^v 



• a topology consistent with a top quark decay: the combination of one 6-tagged jet (b) 
and two untagged jets (j) with the highest p^ b must satisfy mjjb G [120, 240] GeV. 

For the selected events, the transverse mass wit is defined as: 



m T = a/2?W ss (1 - cos A^ imiss ), (6.1) 
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where A</> T ,miss is the azimuthal angle between the r jet and the direction of the missing 
momentum. This discriminating variable is related to the W boson mass in the W — >■ rv 
background case and to the H + mass for the signal hypothesis. 

6.2 Data-driven estimation of the multi-jet background 

The multi-jet background is estimated by fitting its E™ 1SS shape (and the E™ 1SS shape of 
other backgrounds) to data. In order to study this shape in a data-driven way, a control 
region is defined where the r identification and 6-tagging requirements are modified, i.e. 
r candidates must pass a loose r identification but fail the tight r identification used in 
the signal selection, and the event is required not to contain any 6-tagged jet. Hence, 
the requirement on rrijjb is also removed. Assuming that the shapes of the E™ 1SS and tot 
distributions are the same in the control and signal regions, the _E™ 1SS shape for the multi-jet 
background is measured in the control region, after subtracting the simulated background 
contributions from other processes. These other processes amount to less than 1% of the 
observed events in the control region. The E!f lss shapes obtained with the r+jets selection 
of Section 6.1 or in the control region are compared just before the E!f lss requirement in the 
baseline selection in Fig. 4(a). The differences between the two distributions are accounted 
for as systematic uncertainties. For the baseline selection, the E!f lss distribution measured 
in the data is then fit using two shapes: the multi-jet model and the sum of other processes 
(dominated by tt and VF+jets), for which the shape and the relative normalisation are 
taken from simulation, as shown in Fig. 4(b). The ratio between the numbers of multi-jet 
background events in the control and signal regions enters the likelihood function for the 
signal estimation (see Section 8) as a nuisance parameter while the shape of the multi-jet 
background is measured in the same region after additionally requiring E™ 1SS > 65 GeV. 




Figure 4. (a) Shape of E™ lss in a control region of the data or using the baseline selection, 
after subtracting the expectation from ti, W+jets, and single top quark processes estimated from 
simulation. The distributions are compared just before the E™ SB requirement in the baseline 
selection of Section 6.1, with the exception that, in the control region, the r selection and the 
6-tagging requirements are modified, see text, (b) Fit of the E™ lss template to data, in the signal 
region. Only statistical uncertainties are shown. 
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6.3 Backgrounds with electrons and jets misidentified as r jets 

The methods described in Section 5.3 are used to estimate the probability for electrons 
or jets to be misidentified as r jets. The estimated contribution to the background from 
the jet — > t misidentification after the r+jets selection is given in Table 5. The back- 
grounds arising from the jet — > r misidentification are not expected to be well modelled in 
simulation, which is why they are estimated using data-driven methods. 



Sample 


Data-driven method [events] 


Simulation [events] 


tt 


33 ± 1 


37 ± 1 


W+jets 


2.5 ±0.1 


3.9 ± 1.5 


Single top quark 


1.3 ±0.1 


2.0 ±0.3 



Table 5. Application of the misidentification probability obtained from a control region in the data 
enriched in VF+jets events, for the r+jets channel. The predictions of the background contributions 
based on data-driven misidentification probabilities and on simulation are given, with statistical 
uncertainties only. In both cases, all top quarks decay via t — > bW. 



6.4 Data-driven estimation of backgrounds with correctly reconstructed r jets 

An embedding method [54] is used to estimate the backgrounds that contain correctly 
reconstructed r jets. The method consists of selecting a control sample of it-like /i±jets 
events and replacing the detector signature of the muon by a simulated hadronic r decay. 
These new hybrid events are then used for the background prediction. In order to select 
this control sample from the data, the following event selection is applied: 

• event triggered by a single-muon trigger with a pt threshold of 18 GeV; 

• exactly one isolated muon withpx > 25 GeV, no isolated electron with > 20 GeV; 

• at least four jets with pr > 20 GeV, at least one of which is 6-tagged; 

• E™ ss > 35 GeV. 

This selection is looser than the selection defined in Section 6.1 in order not to bias 
the control sample. The impurity from the background with muons produced in r decays 
and non-isolated muons (dominantly bb and cc events) is about 10%. However, this contri- 
bution is greatly reduced as these events are much less likely to pass the r±jets selection, 
in particular the p^ requirement. 

The shape of the tot distribution for the backgrounds with true r jets is taken from 
the distribution obtained with the embedded events, after having applied the r±jets event 
selection. The normalisation is then derived from the number of embedded events: 

gT+_B™ lss — trigger 

N T = Vcmbedded • (1 - c T ^) — /f _ IDitriggcr ■ B( T ~> hadrons + v), (6.2) 
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where N T is the estimated number of events with correctly reconstructed r jets, iV em b e dded is 
the number of embedded events in the signal region, c T ^^ is the fraction of events in which 
the selected muon is a decay product of a r lepton (taken from simulation), e r +- E T lsS-tri gg er 
is the t + E™ 1SS trigger efficiency (as a function of p\ and E™ 1SS , derived from data), 
e ^-iD,tngger j g ^ e muon trigger and identification efficiency (as a function of px and r], 
derived from data) and B{r — > hadrons + v) is the branching ratio of the r lepton decays 
involving hadrons. The hit distribution for correctly reconstructed r jets, as predicted by 
the embedding method, is shown in Fig. 5 and compared to simulation. 




Figure 5. Comparison of the wit distribution for correctly reconstructed r jets, predicted by the 
embedding method and simulation. Combined statistical and systematic uncertainties (as described 
in Section 7) are shown. 



6.5 Event yields and distribution after the selection cuts 

Table 6 shows the expected number of background events for the SM-only hypothesis and 
the observation in the data. The total number of predicted events (signal+background) 
in the presence of a 130 GeV charged Higgs boson with B(t — > bH + ) = 5% is also shown. 
The number of events with a correctly reconstructed r jet is derived from the number of 
embedded events and does not depend on the cross section of the tt — > bbW + W~ process. 
On the other hand, the r+jets analysis relies on the theoretical inclusive tt production 
cross section a t i = 167^Jg pb [26] for the estimation of the background with electrons or 
jets misidentified as r jets. In the presence of a charged Higgs boson in the top quark 
decays, with a branching ratio B{t — > bH + ), the contributions of tt — > bbW + W~ events in 
these backgrounds are scaled according to this branching ratio. The data are found to be 
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consistent with the estimation of the SM background. The ijit distribution for the r+jets 
channel, after all selection cuts are applied, is shown in Fig. 6. 



Sample 


Event yield (r+jets) 


True r (embedding method) 


210 ± 10 ± 44 


Misidentified jet — > r 


36 ± 6 ±10 


Misidentified e — > r 


3± 1± 1 


Multi-jet processes 


74 ± 3 ±47 


All SM backgrounds 


330 ± 12 ± 65 


Data 


355 


t -> bH + (130 GeV) 


220 ± 6 ±56 


Signal+background 


540 ± 13 ± 85 



Table 6. Expected event yields after all selection cuts in the r+jets channel and comparison with 
4.6 fb -1 of data. The numbers in the last two rows, obtained for a hypothetical H + signal with 
m H + = 130 GeV, are obtained with B(t —> bH + ) = 5%. The rows for the backgrounds with 
misidentified objects assume B(t — > bW) = 100%. Both statistical and systematic uncertainties are 
shown, in this order. 




Figure 6. Distribution of mi after all selection cuts in the r+jets channel. The dashed line 
corresponds to the SM-only hypothesis and the hatched area around it shows the total uncertainty 
for the SM backgrounds. The solid line shows the predicted contribution of signal+background 
in the presence of a charged Higgs boson with m H + = 130 GeV, assuming B(t — > bH + ) = 5% 
and B(H + — > tv) = 100%. The contributions of tt — s- bb\V + W~ events in the backgrounds with 
misidentified objects are scaled down accordingly. 
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7 Systematic uncertainties 



7.1 Systematic uncertainties arising from the detector simulation 

Systematic uncertainties arising from the simulation of pile-up and object reconstruction 
are considered. The latter arise from the simulation of the trigger, from the reconstruction 
and identification efficiencies, as well as from the energy/momentum scale and resolution 
for the objects described in Section 3. To assess the impact of most sources of systematic 
uncertainty, the selection cuts for each analysis are re-applied after shifting a particular 
parameter by its ±1 standard deviation uncertainty. The systematic uncertainties related 
to the electrons and muons are discussed in, respectively, Ref. [42] and Refs. [43, 55]. For the 
jets, see Ref. [48] and, in particular, Ref. [50] for the 6-tagging calibration. The systematic 
uncertainties related to r jets are discussed in Ref. [51]. Finally, for the reconstruction of 
i?™ lss , see Ref. [52]. All studies of systematic uncertainties have been updated with the 
full dataset collected in 2011. The dominant instrumental systematic uncertainties arise 
from the jet energy resolution (10-30%, depending on arid 77), the jet energy scale (up 
to 14%, depending on and rj, to which a pile-up term of 2-7% and a b jet term of 2.5% 
are added in quadrature), as well as the 6-tagging efficiency (5-17%, depending on and 
77) and misidentification probability (12-21%, depending on pt and 77). In comparison, 
the systematic uncertainties arising from the reconstruction and identification of electrons 
and muons are small. All instrumental systematic uncertainties are also propagated to the 
reconstructed £™ 1SS . 

7.2 Systematic uncertainties arising from the generation of tt events 

In order to estimate the systematic uncertainties arising from the tt generation and the 
parton shower model, the acceptance is computed for tt events produced with MC@NLO 
interfaced to HERWIG/ JIMMY and POWHEG [56] interfaced to PYTHIA. For the signal 
samples, which are generated with PYTHIA (i.e. without higher-order corrections), no 
alternative generator is available. Instead, the systematic uncertainty for the signal samples 
is set to the relative difference in acceptance between tt events generated with MC@NLO 
interfaced to HERWIG/ JIMMY and with AcerMC, which is also a leading-order generator, 
interfaced to PYTHIA. The systematic uncertainties arising from initial and final state 
radiation are computed using tt samples generated with AcerMC interfaced to PYTHIA, 
where initial and final state radiation parameters are set to a range of values not excluded 
by the experimental data [57]. The largest relative differences with respect to the reference 
sample after full event selections are used as systematic uncertainties. The systematic 
uncertainties arising from the modelling of the tt event generation and the parton shower, 
as well as initial and final state radiation, are summarised in Table 7 for each analysis. 

7.3 Systematic uncertainties arising from data-driven background estimates 

The systematic uncertainties arising from the data-driven methods used to estimate the 
various backgrounds are summarised in Table 8, for each of the three channels considered 
in the analysis. 
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Source of uncertainty 


Normalisation uncertainty 


lepton+jets: 




VJCllCldljUl £llll_l JJCll LUll ollUWCl l uUf r ±1 . OlgilCll ICglvJIll 

Generator and parton shower (bbW~^~ W~ , signal region) 
Generator and parton shower (bbWH + , control region) 
Generator and parton shower (bbW + W~ , control region) 
Initial and final state radiation (signal region) 
Initial and final state radiation (control region) 


10% 

8% 
7% 
6% 
8% 
13% 


r+lepton: 




Generator and parton shower (bbWH + ) 
Generator and parton shower (bbW + W~) 
Initial and final state radiation 


2% 
5% 
13% 


r+jets: 




Generator and parton shower (bbWH + ) 
Generator and parton shower (bbW + W~) 
Initial and final state radiation 

Table 7. Systematic uncertainties arising from the modelling of tt 
events and the parton shower, as well as from initial and final stat 

For backgrounds with misidentified leptons, discussed 


5% 
5% 
19% 

-> bbW+W- and it -> bbWH+ 
2 radiation. 

m Sections 4.2 and 5.2, the 



main systematic uncertainties arise from the simulated samples used for subtracting true 
leptons in the determination of the misidentification probabilities. These are sensitive to 
the instrumental systematic uncertainties and to the sample dependence (misidentification 
probabilities are calculated in a control region dominated by gluon-initiated events, but 
later used in a data sample with a higher fraction of quark-initiated events). 

The dominant systematic uncertainties in the estimation of the multi-jet background in 
the r+jets channel, described in Section 6.2, are the statistical uncertainty of the fit due to 
the limited size of the data control sample and uncertainties due to potential differences of 
the E™ 1SS shape in the signal and control regions. The dominant systematic uncertainties in 
estimating the contribution of events with electrons misidentified as r jets in Sections 5.3 
and 6.3 arise from the subtraction of the multi-jet and electroweak backgrounds in the 
control region enriched with Z — > ee events and from potential correlations in the selections 
of the tag and probe electrons. For the estimation of backgrounds with jets misidentified 
as hadronically decaying r leptons, also discussed in Sections 5.3 and 6.3, the dominant 
systematic uncertainties on the misidentification probability are the statistical uncertainty 
due to the limited control sample size and uncertainties due to the difference of the jet 
composition (gluon or quark-initiated) in the control and signal regions, which is estimated 
using simulation. Other uncertainties come from the impurities arising from multi-jet 
background events and from true hadronic r decays in the control sample. The systematic 
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Source of uncertainty 


Normalisation uncertainty 


Shape uncertainty 


lepton-)-jets! lepton misidentification 




Choice of control region 


6% 




Z mass window 


4% 




Jet energy scale 


16% 




Jet energy resolution 


7% 




Sample composition 


31% 




T-t-lpntofi ' let — ^ t Tnisifipntification 




Statistics in control region 


2% 




Tpt pom ti pitii "1" l nn 


11% 




Oni ppt-vpI afpn t;"\rt;'l"P'm a 1 1 p*3 
V-zujcv^ij i cicxtC'V-L oy o tciiicxijiv^o 


23% 


3% 


t 1 1 p "fit on ' p V t m i qi H p-nti ti Pti ti p»n 

I ~1CU Umil. O T / llllolLlCllLlllLjCLLlCll 




\/l ieiriPi'vl~'i'riP£i'l~mn T~>T"r^r^cir^ili'l"\7' 
IvllolU.cllUlllLclljlL'Il UI UUdUlll Ij V 


20% 




T-r-lpnton " Ipnton m i si Hpntifi ration 




Choice of control region 


4% 




Z mass window 


5% 




Jet energy scale 


14% 




Jet energy resolution 


4% 




Sample composition 


39% 








H ,TT1 P» P/H H 1 Tl fl" novo m P"l~PT*c; 

11/IIlUcU.U.lllg pal cLillcltJI t> 


6% 


3% 


iviuoii looiaLioii 


7% 


2% 


PavainptpT''; in Tioi'TTifilitiatinn 

J. CLL (X111C ICIO 111 1HJ1 llldllOdlilUll 


16% 




t lHpntiripnf inn 

/ H_11^111j11H_/CIIj1U11 


5% 




t pn py o~\t era 1 p 
/ c-nc-i gy ov^cxic 


6% 


1% 


t— 1— iptci" ipt — V -T in mi n pnt l fi pa t inn 

/ ^^ICIjO. ICl T / llllolvJ-CllulllvjCHjlwll 




St^tmtips m pnrit"rnl rpcnnn 

kJ LicliLilD IjIvjO 111 vjVJIIuI \JL 1 CC1U11 


2% 




Tpt poTTrnnsitinn 

U vU *^*_/lll LHJ11 


12% 




Mi i vi l ti rnntrn vpcnpin 

1 U.1 1 111 i^UULl Ul I CglUll 


6% 


1% 


Object-related systematics 


21% 


2% 


r+jets: e — > r misidentification 




Misidentification probability 


22% 




r+jets: multi-jet estimate 




Fit-related uncertainties 


32% 




E'n^-shape - m con t ro l region 


16% 





Table 8. Dominant systematic uncertainties on the data-driven estimates. The shape uncertainty 
given is the relative shift of the mean value of the final discriminant distribution. A "-" in the 
second column indicates negligible shape uncertainties. 
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uncertainties affecting the estimation of the background from correctly reconstructed r 
jets in the r+jets channel, discussed in Section 6.4, consist of the potential bias introduced 
by the embedding method itself, uncertainties from the trigger efficiency measurement, 
uncertainties associated to simulated r jets (r energy scale and identification efficiency) 
and uncertainties on the normalisation, which are dominated by the statistical uncertainty 
of the selected control sample and the r + E™ 1SS trigger efficiency uncertainties. 

8 Results 

In order to test the compatibility of the data with background-only and signal+background 
hypotheses, a profile likelihood ratio [58] is used with (lepton+jets), E™ 1SS (r+lepton) 
and rriT (r+jets) as the discriminating variables. The statistical analysis is based on a 
binned likelihood function for these distributions. The systematic uncertainties in shape 
and normalisation are incorporated via nuisance parameters, and the one-sided profile 
likelihood ratio, q^, is used as a test statistic. No significant deviation from the SM 
prediction is observed in any of the investigated final states in 4.6 fb _1 of data. Exclusion 
limits are set on the branching fraction B(t — >■ bH + ) and, in the context of the m™ ax scenario 
of the MSSM, on tan /3, by rejecting the signal hypothesis at the 95% confidence level (CL) 
using the CL S procedure [59]. These limits are based on the asymptotic distribution of 
the test statistic [58]. The combined limit is derived from the product of the individual 
likelihoods, and systematic uncertainties are treated as correlated where appropriate. The 
exclusion limits for the individual channels, as well as the combined limit, are shown in 
Fig. 7 in terms of B(t — > bH + ) with the assumption B(H + —?■ tu) = 100%. In Fig. 8, the 
combined limit on B{t — > bH + ) x B(H + —> tu) is interpreted in the context of the m™ ax 
scenario of the MSSM. The following relative theoretical uncertainties on B(t — > bH + ) are 
considered [60, 61]: 5% for one-loop electroweak corrections missing from the calculations, 
2% for missing two-loop QCD corrections, and about 1% (depending on tan/3) for Aft- 
induced uncertainties, where A5 is a correction factor to the running b quark mass [62]. 
These uncertainties are added linearly, as recommended by the LHC Higgs cross section 
working group [61]. 

9 Conclusions 

Charged Higgs bosons have been searched for in ti events, in the decay mode t — > bH + fol- 
lowed by H + tv. For this purpose, a total of 4.6 fb _1 of pp collision data at yfs = 7 TeV, 
recorded in 2011 with the ATLAS experiment, is used. Three final states are considered, 
which are characterised by the presence of a leptonic or hadronic r decay, E™ 1SS , b jets, and 
a leptonically or hadronically decaying W boson. Data-driven methods and simulation are 
employed to estimate the number of background events. The observed data are found to 
be in agreement with the SM predictions. Assuming B{H + — > tu) = 100%, upper limits 
at the 95% confidence level have been set on the branching ratio B{t — > bH + ) between 
5% (m,ff+ = 90 GeV) and 1% (m H + = 160 GeV). This result constitutes a significant 
improvement compared to existing limits provided by the Tevatron experiments [14, 15] 
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Figure 7. Expected and observed 95% CL exclusion limits on B(t bH + ) for charged Higgs 
boson production from top quark decays as a function of m H +, assuming B(H + — > tv) = 100%. 
Shown are the results for: (a) lepton+jets channel; (b) r+lepton channel; (c) r+jets channel; (d) 
combination. 



over the whole investigated mass range, but in particular for m#+ close to the top quark 
mass. Interpreted in the context of the m™ ax scenario of the MSSM, tan (3 above 12-26, as 
well as between 1 and 2-6, can be excluded in the mass range 90 GeV < mjj+ < 150 GeV. 
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